Movers and shakers: Granular damping in microgravity 
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The response of an oscillating granular damper to an initial perturbation is studied using experiments per- 
formed in microgravity and granular dynamics simulations. High-speed video and image processing techniques 
are used to extract experimental data. An inelastic hard sphere model is developed to perform simulations 
and the results are in excellent agreement with the experiments. The granular damper behaves like a frictional 
damper and a linear decay of the amplitude is observed. This is true even for the simulation model, where 
friction forces are absent. A simple expression is developed which predicts the optimal damping conditions for 
a given amplitude and is independent of the oscillation frequency and particle inelasticities. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

The characteristic property of dynamic granular systems, 
when compared to other many-particle systems, is their abil- 
ity to dissipate mechanical energy through particle collisions. 
While the dissipative properties of vibrated granulate have 
long been investigated (H, 0], recently a large body of liter- 
ature OHS!] has emerged on the mechanics and technical ap- 
plication of this damping mechanism in the form of granular 
dampers. A granular damper is a container partly filled by 
granular particles which may be attached to vibrating machin- 
ery to attenuate the amplitude of the oscillations. In its regime 
of operation, the granular material is in a gaseous state and 
its dynamics is determined primarily by the interparticle colli- 
sions rather than by long-lasting sliding contacts between the 
grains. Static granular dampers (e.g., Refs. [1(3]) which 
exploit the rheology of granular matter and impact dampers 
(e.g., Refs. IUHmI), where only one or few particles are lo- 
cated in a cavity and dissipate energy in collisions with the 
walls of the container, are not considered here. 

Granular dampers have a number of properties which are 
desirable in a wide range of technical applications: Unlike tra- 
ditional dampers, granular dampers do not require an anchor 
in order to restrict the motion of the system. This is advanta- 
geous for damping in portable equipment and in space appli- 
cations where no fixed anchor is available. Granular dampers 
are extremely simple devices consisting solely of particles en- 
closed in a container or cavity and require very little mainte- 
nance. Granular dampers do not suffer from significant aging 
when compared to the oil and rubber components of tradi- 
tional dampers. Finally, granular dampers can operate over 
a wide range of temperatures without performance degrada- 
tion as the mechanics of the particle-particle and particle-wall 
interactions exhibit only a weak dependence on the tempera- 
ture. Modern technical applications of granular dampers in- 
clude the damping of blade integrated disks (blisks) for com- 
pressors 11511 . structural vibration damping [ 16 -} 1811. noise re- 
duction of bank note processing machines [19] and others. 
Perhaps the most common application is the dead-blow ham- 
mer l2(ill and in other impact damping handles I2H1 . 

The macroscopic damping properties of granular dampers 
under dynamic load is complicated, highly non-linear, and 
there is no straightforward way to optimize their performance 



for a given situation. This has been demonstrated in a number 
of experiments and Molecular Dynamics (MD) simulations, 
including investigations on the attenuation of a fre e sp ring 
or cantilever with an attached granular damper (§L I7L l22l - l24tl . 
The response of an oscillating cantilever with respect to peri- 
odic forcing has also been studied f25T430ll . Even more com- 
plex systems have been investigated, such as the oscillation 
modes of a plate with an abundance of granulate filled cav- 
ities i3ll - [34ll with the aim of noise reduction lfl9tl . For sim- 
ple systems, such as cantilever oscillators, some progress has 
been made. Theoretical models have been developed based on 
phenomenological descriptions of the multiphase gas-particle 
flow of granular matter for attenuating oscillations IB5I1 and 
also for driven steady state oscillations 113611 . 

A granular damper, that is a dynamical system of dissipa- 
tive interacting particles, obviously must be able to dissipate 
energy; however, its general behavior is not clear a priori. 
Properties, such as the dissipation rate, are complex functions 
of the frequency and amplitude of the oscillation, as well as 
the particle properties, the extension and characteristics of the 
container or cavities, and the filling fractions. More work is 
needed in this field to generate experimental results and corre- 
sponding models capable of describing the dynamics of gran- 
ular dampers. 

Saluena et al. yj] have shown that several regimes of energy 
dissipation exist for a granular damper and that the transitions 
between these regimes are determined primarily by the influ- 
ence of gravity. An efficient operation of a granular damper 
can be expected only if the average kinetic energy of the parti- 
cles is much larger than their average potential energy and the 
damper operates in the dynamic or collisional regime. In or- 
der to carefully investigate this regime, the influence of gravity 
should be minimized and experimental investigations should 
be performed under conditions of weightlessness. 

The objective of this paper is to develop an effective model 
for the energy dissipation of a granular damper operating in 
the collisional regime. Our approach is as follows: First, 
experiments in microgravity are performed and the attenua- 
tion of a spring with an attached granular damper for sev- 
eral sets of parameters is obtained (Sec. [H). A model capa- 
ble of reproducing the experimental results is also developed 
and high-precision Discrete Element Method (DEM) simula- 
tions are performed (Sec. [TTTb - The two free parameters of 
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the model (inelasticities) are obtained by adjusting the values 
until the simulation matches the experiment as closely as pos- 
sible for a single experiment (Sec. IIII CI . From the excellent 
agreement of the simulation results for the fitted system and 
for all other experiments, it is concluded that the model un- 
derlying the simulation replicates the system's essential fea- 
tures (Sec. IIIIDI i. Thus, the DEM simulations are an effec- 
tive model for granular damping in the collisional regime. In 
Sec. [IV] a simple equation for the optimal design of a simple 
damper is derived and tested against the results of the DEM 
simulations. Section[V]discusses the observed linear decay of 
the amplitude and compares it to friction-damped oscillators. 
Finally, in Sec.[VT]the conclusions of the paper are outlined. 



II. EXPERIMENTAL SETUP 

Figure[T]is a diagram of the experimental setup. Our gran- 
ular damper comprises of a container of adjustable length 
which is partially filled with granular material. The damper 
is mounted to one end of a spring-steel blade and the opposite 
end is clamped in a solid aluminum base plate. The spring 
blade is described fully in Sec. IIII Cl The rectangular damper 
container is constructed from 5 mm thick transparent poly- 
carbonate plates. The internal dimensions of the container 
are 50 mmx50 mmxL, where the length L (in the direction 
of the oscillation) is adjusted by altering the spacing of the 
end walls. The container's net weight (without granulate) is 
M = 434 g. In this work, four different container lengths of 
L = 40, 65, 85, and 104 mm are used. The damper is loaded 
with 37 precision steel ball-bearings of diameter a = 10 mm 
and mass m = 4.04 g. This number of particles is chosen as 
it packs to form a dense layer, two particles deep, on the end 
walls of the container. 
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FIG. 1: Schematic of the experimental setup in front view (left) and 
side view (right). The curvature of the oscillations is exaggerated for 
the purpose of illustration. 

The motion of the damper and contained granulate is 



recorded using a high speed camera (MotionScope M3™), 
which records at a frame rate of 500 fps and with a spatial res- 
olution of 1024 x 1280 pixels. A 45° mirror is placed at the 
side of the container and allows for the simultaneous observa- 
tion of the granulate from the top and the side. The position 
of the damper and the center of gravity of the particles are 
extracted from the top-view using standard image processing 
techniques. Although the side view facilitates more complex 
methods of reconstruction, it will be shown that the motion of 
the granulate is well described by the center of mass motion. 
All position measurements are made in a 2D plane which is 
parallel to and intersecting the top of the container while it is 
in its equilibrium position. Although this introduces some er- 
ror at the peaks of the oscillations due to the curved damper 
trajectory, this error is negligible. 

At the start of the experiment, the spring blade is deflected 
and held at the initial displacement of An = 107.5 mm using 
an electromagnet. A trigger mechanism begins the experiment 
and starts the camera recording. After a short delay of 1 s, 
the spring blade is released from the electromagnet and the 
oscillations are recorded for 30 seconds. 

To assure conditions of weightlessness, the experiment is 
performed on a modified Airbus A300 aircraft which has been 
retrofitted for performing parabolic flights. The parabolic 
flight provides a suitable microgravity environment (±0.05 g) 
which lasts around 22 s and allows a number of experiments to 
be performed. In the following section, the numerical model 
and simulation techniques are described. 



III. NUMERICAL MODEL AND SIMULATION METHOD 

A sufficiently complex model, capable of reproducing the 
observed experimental behavior, must be found if the system's 
dynamics is to be understood. The model presented here is 
complex enough to yield quantitative agreement with the ex- 
periments and yet simple enough to gain insight into the dy- 
namics of the system. The model for the granulate consists 
of a system of N smooth inelastic hard spheres, each of mass 
to and diameter a. Although inelastic hard spheres are a basic 
model for the granulate they capture all of the essential behav- 
ior of the system: dissipative interactions between hard spher- 
ical particles. Friction forces, which typically play a large role 
in granular systems, are also neglected and the implications of 
this approximation are discussed later. 

To model the oscillating mass and granular damper, the 
hard spheres are shaken in a rectangular box of mass M, 
which is coupled to a linear spring. The box is assumed to 
remain parallel to the axis of the system and oscillate along 
only one axis. By only permitting oscillations in a single di- 
mension, this model neglects the arcing motion of the blade 
spring (see Fig. [U and oscillatory modes induced by the col- 
lisions of the particles with the box; however, these effects 
are expected to be small in comparison to the dynamics of the 
modeled oscillation. With these approximations, the collision- 
free motion of the box can then be modeled using a simple 
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harmonic oscillator 

r Box (t) = n A cos [2 vr w (i + Uift )] + ^ ( 1 ) 



where reox is the current position of the oscillator, fj^J x is its 
equilibrium position, n is the unit vector in the direction of 
the oscillation, and ui is the frequency of the empty damper. 
The amplitude of the oscillation A and the phase shift of the 
oscillator t s hift are dynamical quantities and are altered by 
particle-box interactions. At time t = 0, the plate is at its 
positive maximum (i S Mft = 0) with an initial displacement of 
A = A . 

The methods for performing event-driven simulations using 
smooth hard-spheres and fixed walls are well established (e.g., 
see Ref. 13711 ) and will not be discussed in detail. Here an 
event-driven dynamics simulation package (DYNAMO I38IP is 
used to simulate the dynamics of granular-damped oscillators. 
The only extension to the basic event-driven method concerns 
the detection and execution of events between particles and 
the oscillating boundary walls perpendicular to the oscillation 
direction h, which is discussed in the following subsections. 



A. Detecting Oscillating Wall Interactions 

Event-driven algorithms require an expression to calculate 
if and when a collision (an event) occurs between a particle 
and the bounding walls of the damper. If a collision is detected 
and it is the next event to occur in the system, the system is 
moved to the time of the collision and the event is executed by 
updating the velocities of the colliding particle, and the phase 
shift t B hift and amplitude A of the oscillator. 

To determine the time at which a particle i will collide with 
an oscillating wall, the equations of motion for the particle 
and the oscillating plate must be solved. Essentially, this is a 
search for the shortest positive root At of the function 



/(At) = [n + Atv % - r Box (Ai + t)]-n± 



L-a 



-- 
(2) 



where r*i and v\ are the position and velocity of particle i at 
the current system time t, and At is the time till collision. The 
sign of the term ±(L — u)/2 is used to set which side of the 
oscillating box is tested for collisions. 

To guarantee that no roots are missed, the root finding tech- 
nique of Frenkel and Maguire [39] is used. This root finding 
routine requires a fixed interval to search for possible roots. 
The upper bound on the interval to search is determined from 
the time the freely moving particle takes to cross the extrema 
of the tested-wall's oscillation 
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n(n • Vi) ([L - a] /2 + A) - h • (r, ; - r^) 



(3) 



where sgn(ir) is the sign function. The lower bound is typi- 
cally the current system time (Ai m ; n = 0); however, if the last 



event to occur was a collision between this particle and a oscil- 
lating wall, the lower bound is increased to avoid re-detecting 
the same root. The lower bound is then set to 



Ai m ; n — 



2/(0) 



fn 



(4) 



where / max = A ui 2 is the maximum absolute second deriva- 
tive of Eq. (f2|i. The root finding technique used to search 
for suitable roots of Eq. iterates towards a root from the 
boundaries of the interval by approximating the function at 
each iteration with a parabola. The equation of the parabola is 
generated using the derivatives of Eq. (O and its smallest root 
provides the next iteration point. The iterations are halted on 
the nth iteration once the following criterion is met 



\At n - Atn_i| < [Z/ . g] x 1(T 12 

2 /max 



(5) 



where / max = \v{ ■ n| +2 7ru; A is the maximum absolute first 
derivative of Eq. (0. Unlike the hard line system of Frenkel 
and Maguire H39I1 . all roots of Eq. (O are acceptable and only 
the earliest root must be found. This completes the description 
of the collision detection and root finding technique. 



B. Executing Particle-Oscillating Wall Collisions 

The final part of the simulation algorithm concerns the ex- 
ecution of oscillating boundary wall collisions. The conser- 
vation of momentum and the assumption of a constant coeffi- 
cient of inelasticity leads to 



Apt = -Ap waU = -fi(l + e 



pw) (j 



h (6) 



where Api and Ap wa n are the momentum change of the 
colliding particle i and oscillating wall respectively, p, = 
m M/{m + M) is the relative mass, and e pw is the coefficient 
of inelasticity for particle-wall collisions. During a collision, 
the phase, t s hift, and amplitude, A, of the oscillating wall are 
altered under the constraints of conserving momentum and the 
current box position. This results in the following expressions 
for the post collision state of the oscillating plate 
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(8) 



where the primes denote post-collision values. Care must be 
taken at this point in the calculation to ensure that the magni- 
tude of t and t s hift do not affect the precision of the calcula- 
tions. Care must be taken also to retain the correct quadrant 
of the calculated angle when using the arctan function. 

A difficulty with the event-driven simulation method arises 
from its inability to simulate events with finite durations. 
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When the oscillating wall is accelerating, a particle can re- 
peatedly collide with the plate until its relative velocity and 
separation are numerically zero. Physically, the particle sticks 
to the wall and is pushed until the plate enters the deceleration 
phase of its oscillation, or interacts with another particle. To 
prevent this unresolvable situation from occurring within the 
event-driven simulation, the interactions between the oscillat- 
ing wall and a particle are turned elastic when 
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< 0.04 



(9) 



The pushing of the particle is then transformed into a sequence 
of small hops which, as in the physical pushed case, do not 
dissipate energy. As this expression is linear in the current 
displacement A, the long time behavior of the system is still 
recovered (A — > as t — > oo). This elastic approximation 
is small when the plate motion dominates the dynamics of 
the system and the results appear to be unaffected if smaller 
values for Eq. (O are used. 



C. Parameters of the Simulation 

The simulations are initialized with all particles arranged in 
a regular lattice (FCC), with initial velocities assigned from 
a Gaussian and a total particle energy less than < 0.002% of 
the initial spring energy. The particles are packed in a loose 
layer on the wall at the initial extrema of the oscillation. The 
particles in the experiment are also typically arranged this way 
due to the influence of gravity before the microgravity phase 
of the experiment. 

The simulations require several inputs and these parame- 
ters are reported in Table U All parameters, with the excep- 
tion of the box frequency cj and inelasticities e pp and e pw , 
are directly obtained from the experimental setup described in 
Sec. Q]] The three remaining parameters must be calculated 
from material parameters or obtained through experimental 
results. 



TABLE I: Model parameters for the event-driven simulations. 



a (mm) 


m (g) 


N 


Ao (mm) 


u^s" 1 ) 


M(g) 


Epp 




10 


4.04 


37 


107.5 


1.23 


434 


0.75 


0.76 



The frequency of the unloaded damper u may be estimated 
using the simple harmonic oscillator model. The spring con- 
stant of the spring-blade may be calculated using the Euler- 
Bernoulli beam equation 



k = 



Ewh 3 
4/ 3 



= 0.0254 N mm" 



(10) 



where E — 2.06 x 10 5 N mm" 2 is the elastic modulus of the 
spring steel, w = 30 mm is the spring width, h — 1.5 mm 
is the spring thickness, and I = 590 mm is the spring length. 
If the system behaves as a simple harmonic oscillator and the 



mass of the spring is ignored the frequency may be estimated 
using 



1 k 
^ ~ ^Vtt ~ l-217s" 
2tt V M 



(11) 



The frequency of the loaded granular damper (w sys tcm) is 
lower than that of the empty damper (lj) due to the added mass 
and the interactions of the granulate. In the simple harmonic 
oscillator model, the additional mass of the granulate alters 
the frequency of the oscillations by 



^loaded — W 



M 



M + Nm 



(12) 



In the limit that the granulate is tightly packed in the gran- 
ular damper, the frequency of the system should limit to the 
simple harmonic oscillator frequency w sys t e m — >• wi oa dcd- In 
the limit of a large box, the granulate will completely decou- 
ple from the oscillator and w sys t e m — > w. Remarkably, the 
frequency of the experimental oscillators, obtained through 
averaging the peak and center point frequencies, is consis- 
tent for all box lengths at approximately w S ystcm ~ 1.05 s — 1 
with a standard deviation of ±0.01 s _1 . If its assumed that 
^system ~ ^loaded for small box lengths, Eq. dT2b estimates an 
unloaded frequency of U) rj 1.22 ± 0.1 s _1 for the experimen- 
tal system. This agreement with the beam equation is promis- 
ing and suggests that, although the granulate is periodically 
decoupled from the oscillator, the deviation from Eq. ( TTZb is 
still small for the experimental box lengths studied here. For 
the simulations, a slightly higher frequency of to « 1.23 s _1 
is used which is within the standard deviation of the experi- 
mental values and yields an excellent fit to the experimental 
data. 

describ- 



Finally, the coefficients of restitution e pp and 



pw 



ing the inelastic collisions between particles and between a 
particle and the wall must be determined. These model pa- 
rameters are obtained by fitting simulation results to the exper- 
imental data for the smallest box length (L = 40 mm, Fig. [2]). 
As best fits the following results are obtained 



0.75 



0.76 



(13) 



The value for the particle-wall coefficient of restitution is in 
close agreement with published results reported for a 9.35mm 
steel ball-bearing impacting a clamped acrylic plate i40ll ; 
however, the particle-particle inelastic coefficient is signifi- 
cantlylower than expected. Performing an automated drop 
test H41TI of the granulate on to a silicon carbide plate yields 
an elasticity of e ~ 0.95. Due to the high rigidity of the base 
plate, this value should be close to the experimental value for 
particle-particle interactions. The fitted particle-particle in- 
elasticity £ pp may be unexpectedly lower than the drop test 
results due to missing dissipation mechanisms in the model 
(e.g., granulate friction). Despite this, the agreement of the 
simulation and experimental results (see Sec. lIIlDl i shows that 
this is still an effective model for the system. 

It should be noted that the optimization/fitting of the in- 
elasticities £ pp and e pw is performed exclusively for the box 
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width of L — 40 mm. For all other simulations reported here, 
the optimal coefficients of restitution are used without further 
fitting. 

D. Validation of the Numerical Method 

The simulation and experimental results are compared in 
this section to validate the model. Figure [2] presents the box 
position xbox and granulate center of mass cccom as a func- 
tion of time for a box length of L = 40 mm. Two experimen- 
tal measurements are reported and both are in close agreement 
with the simulation results. The experimental and simulation 
results display a high degree of repeatability and single real- 
izations are representative of the averaged values. This is due 
to the uniqueness of the initial state, with the spring held in a 
deflected state and the particles resting in a regular, repeatable 
layer on the outer wall due to the influence of gravity before 
the microgravity phase. However, the experimental results be- 
gin to fluctuate towards the end of the microgravity phase due 
to disturbances in the flight. 




FIG. 2: (Color online) A comparison of simulation results, experi- 
mental data and Eq. dl7b for (a) the box position ibox and (b) the 
granulate center of mass icom as a function of time for a box length 
of L — 40 mm. The simulation data is fitted to the experimental data 
through the inelasticities e pp = 0.75 and e pw = 0.76. 



The numerical result for the box position xbox as a function 
of time is in excellent quantitative agreement with the experi- 
mental data. For the position of the center of mass, xqom the 
agreement is also very good albeit not as close as for XBox{t), 
with some over-estimations near the peaks of the oscillations. 
The error could arise from the experimental method due to the 
top-down view of the simulation and 2D image reconstruction 
used. The area of the visible particles are identified and the 
centroid location is taken to be the center of mass. Due to 
the end walls and slight arcing motion of the box (see Fig. |T|i 
the reconstructed center of mass is slightly biased towards the 
center of the box. 

The agreement between the simulation and experiment for 
the frequency of the damped oscillator is excellent and con- 
firms the accuracy of the fundamental frequency ui; however, 



the excellent agreement in the amplitudes between experimen- 
tal data and simulations for L = 40 mm is perhaps not too 
surprising since this experimental data set is used to deter- 
mine the coefficients of restitution, e pp and e pw . The model 
parameters are now fixed and the numerical result for several 
different box widths are compared with the corresponding ex- 
perimental data (see Figs. [3][5]). 




FIG. 3: (Color online) A comparison of simulation results, experi- 
mental data and Eq. l !17b for (a) the box position ibox and (b) the 
granulate center of mass icom as a function of time for a box length 
of L — 40 mm. Line types are described in Fig. [2] The simulation 
data is not fitted to the experiment and the parameters of Fig. [2] are 
used. 




FIG. 4: (Color online) The same comparison as Fig. [3] but for a box 
length of L = 85 mm. 

In general, the simulation results are in excellent agreement 
with the models predictions. This implies that the approxima- 
tions of the model (ID oscillations, no air resistance, ideal 
spring) are small and have little effect on the dynamics of the 
granular damper. Some of these approximations may already 
be compensated for in the fitting of the coefficients of resti- 
tution, but they appear to be well behaved with the changes 
in box length. In the simulation, rotational degrees of free- 
dom are neglected by eliminating friction between the parti- 
cles and the particles and the container walls. In contrast to 
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FIG. 5: (Color online) The same comparison as Fig. [3] but for a box 
length of L — 104 mm. 



made to estimate the optimal box length using a simple model 
for the dynamics. 
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FIG. 6: (Color online) Simulation values for the cumulative energy 
loss through interactions with the side walls (P-S), end walls (P-E) 
and particle-particle interactions (P-P) for a box length of L = 40 
mm. 



vibrated granular dampers in gravity, where the energy dissi- 
pation due to friction is of the same scale as energy dissipation 
by impact [42], in microgravity friction seems to be less im- 
portant or easily characterized into the restitution coefficient 
£ pp . Overall, the fitting of the inelasticities appears to be ef- 
fective at capturing the behavior of the system and no further 
parameters or extensions of the simulation model are required. 

The most striking feature of the curves in Figs. [2]|5]is the 
linear decay of the peak amplitude of the oscillation with time. 
A detailed discussion of this property is postponed to Sec.lVl 
and optimal dampers are discussed in the following section. 



IV. OPTIMAL DAMPERS 

There is a significant dependence of the damping efficiency 
on the container length, as is seen in Figs. [2]-|5] The number 
of cycles before the oscillations are sufficiently damped varies 
from 13 to 4 as the box length is increased. By examining the 
energy transfer mechanisms within the granular damper, an 
expression for optimizing the dampers design may be found. 

Figure |6]plots the cumulative energy lost through the three 
classes of collisions in the simulation system. It should be 
noted that Fig. [6] is only valid for the fitted inelasticities and 
will therefore differ from the true experimental values. Never- 
theless, the results should agree qualitatively and allow some 
insight into the experimental system. The sides of the box ap- 
pear to be unimportant in this design of a damper and may 
present an opportunity for optimization by utilizing alterna- 
tive shaker geometries (e.g., an hourglass design). Not only 
are the particle-end wall collisions the sole mechanism for 
the transferal of oscillation energy from the oscillator into the 
contained granulate, but simulation results estimate that these 
collisions are also a significant dissipation mechanism for the 
damper. The end wall interactions both transfer and dissipate 
the maximum energy when the relative velocity of the oscilla- 
tor end walls and granulate are maximized. Therefore, maxi- 
mizing this relative speed should optimize the performance of 
the granular damper. In the following subsection, an attempt is 



A. Theoretical Predictions 

Attempting to optimize the system by modeling the gran- 
ulate as a single particle or some other simplified description 
is difficult [43] due to the lack of an analytical solution to 
Eq. (O. To estimate the optimal damping conditions, only 
two plausible assumptions for the motion of the granulate in 
the box are required: a) the granulate will be "collected" and 
form a packed layer on the approaching oscillating wall dur- 
ing the initial inward stroke (when the oscillator accelerates 
towards the center point), and b) the center of mass velocity 
of the granulate at the end of the inward stroke is, on average, 
the maximum oscillator velocity. The time from the center of 
the stroke at which the granulate would hit the peak displace- 
ment of the oscillator is then given by 



A- 



;,pcak 



^laycr 



2 7T ^loaded A 



(14) 



where <7i ayer = 20 mm is the thickness of the layer of gran- 
ulate when it is packed on the surface of the oscillating wall. 
It should be noted that Eq. ( fT4l decreases in time, as A de- 
creases on average due to interactions with the granulate. If 
for any integer n the peak collision time lies in the range 
n < wig, peak < n + 1 /a, the granulate will hit the oscillat- 
ing wall on the outward phase of its stroke. All experimental 
box with the exception of the largest system (L = 104 mm) 
are within this regime. It is expected that improved damping 
occurs if the granulate hits on the inward stroke as the relative 
velocity is maximized. 

The granulate travels the length of the box in 



L 



^ layer 



2 7T Wi oa( jed A 



(15) 



If n+ 1/4 > u) tg,p C ak > n + 1 and n < 10 i gi Box < n + 1/2, the 
granulate will collide on the inward phase of the stroke. The 
largest system, where to £ gi Box ~ 0.15, collides after the turn- 
ing point of the oscillator; however, the dissipation is maxi- 
mized when co ig.Box ~ V 2 - At this point, the relative velocity 
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between the granulate and oscillating box is also maximized. 
For uj tg.Box > V 2 the pl ate is either decelerating or multiple 
cycles of the oscillation occur without the granulate colliding. 

The damping of the oscillator from the initial state can be 
optimized, independently of the inelastic coefficients, by al- 
tering either L, ui, or A such that ui g ,B x ~ V 2 - Efficiency 
will be lost and recovered as A changes over time, but if the 
granulate is relatively inelastic this will occur after most of the 
energy is dissipated or transferred in the first cycle. 

Setting wt g Box = V 2 m Eq. ( fTSl l and using Eq. ( fI2t . the 
optimal box length L opt may be estimated for a given initial 
amplitude Aq using 



L opt = tt A 0] J — — + a laycr (16) 

This expression is remarkable in that it is independent of the 
oscillation frequency. This may be understood from dimen- 
sional analysis as, due to the negligible initial kinetic energy, 
the model has only one time scale. As such, the solutions to 
the model must scale trivially in the frequency of the oscilla- 
tions. In the following subsection, the results of Eq. ( TToT l and 
its assumptions are checked against simulation results. 

B. Numerical Test 

The validity of the basic assumptions made in Sec. IIV Al 
and the result, Eq. ( TToT l. are now tested using the results of 
the DEM simulations. Using Eq. ( TToT l to predict the optimal 
box length for the damping of the experimental system yields 
a value of L opt = 311 mm. The results of a simulation at this 
box length are presented in Fig. [7] A square step in the gran- 
ulate center of mass velocity is visible at the peak of the box 
velocity as the granulate decouples from the oscillator. The 
assumption of an equal box and granulate velocity at the mid- 
point of the stroke (at peak velocity) appears to hold. Visual 
inspection [44] confirms the granulate is collected in a layer 
on the approaching oscillating wall. The re-collision of the 
granulate also appears to occur close to the peak of the box 
velocity, maximizing the relative velocity, energy dissipation, 
and energy transfer in this first collision. The largest oscilla- 
tions are effectively damped within one second; however, the 
oscillator is now susceptible to smaller amplitude oscillations 
which appear to decay very slowly. The optimal approach 
would be to couple two or more dampers to damp a wider 
range of amplitudes within short timescales. This idea has al- 
ready been pursued for impact dampers (e.g., Ref. [45]) which 
are related to granular dampers except that in the container or 
cavity there is only a single particle. 

To test the predictions of Eq. ( TToT ) for the optimal damping 
length L, a suitable metric must be defined to compare various 
box lengths. Figure[8]compares the time an oscillator takes to 
dissipate a certain fraction of the initial energy as a function of 
the box length. Despite the continuing low-amplitude oscilla- 
tions of the damper at L = 311 mm (see Fig. O, the damper 
effectively eliminates 95% of the initial energy in well under 
two oscillations. Equation (TToT l appears to yield an excellent 
estimate for the global optimal box length, avoiding both the 




t(s) 



FIG. 7: (Color online) A comparison of simulation results, exper- 
imental data and Eq. U7b for the box and granulate (a) position, 
and (b) velocity as a function of time for the optimal box length of 
L — 311 mm, as predicted by Eq. d!6l l. 



highly inefficient zones towards the edges of the graph. Pre- 
vious work on forced granular dampers (e.g., see Fig. 7 in 
Ref. l25ll ) also yields performance curves with the same gen- 
eral U-shape as Fig. [8] An alternative metric for comparing 
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FIG. 8: Simulation results for the time t to dissipate a percentage of 
the initial energy, versus the box length L. The vertical dashed line 
indicates the optimal box length as predicted by Eq. d 1 6b - 

the efficiency of granular dampers is now defined through a 
phenomenological model for the damping behavior. 

V. PHENOMENOLOGICAL MODEL OF GRANULAR 
DAMPERS 

Figures|3f5]reveal a linear decay of the amplitude of the os- 
cillation with time, and thus the energy of the system decays 
quadratically in time. This is highlighted in Fig.[9j where the 
time dependence of the square root of the total energy of the 
damper is plotted. This result is surprising considering the ap- 
proximations of the previous section: the oscillator appears to 
have a constant frequency for a given box length, and the os- 
cillator collects the granulate on a wall and then collides the 
granulate in each half period. For the amplitude decay to be 
linear, the energy dissipated in each of these "collisions" of the 
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granulate must then be proportional to the amplitude A. How- 
ever, if an inelastic particle is given a velocity proportional to 
the maximum plate velocity (2 tt lo A), it will dissipate energy 
proportional to the square of the plate amplitude (A 2 ) for a 
given number of collisions. 
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FIG. 9: (Color online) Total energies of the oscillators as obtained 
from numerical simulations. 



The result is also surprising as the more common viscous 
dampers yield an exponential decay of the amplitude; how- 
ever, the only simple damped-oscillator which displays a lin- 
ear decay in the amplitude is one damped by a constant magni- 
tude force 1 4611 (e.g., a friction-damped oscillator). The equa- 
tion of motion for such an oscillator is 



M tQt x = —kx - n M to t sgn (i) 



(17) 



where x is the oscillator position, M to t = M + N m is the 
total oscillating mass and « M to t is the magnitude of the con- 
stant frictional force. No simple analytical solution exists to 
this equation although piecewise solutions may be found 14611 . 
This model does not appear to be appropriate for the gran- 
ular damper due to the intermittent nature of the damping 
force. For example, the steps in the damper energy (Fig. |9]l 
arise from the decoupling of the damping granulate and oscil- 
lator at the midpoint of the stroke (see Fig. [7]) and during these 
steps the oscillator experiences no damping force. However, 
this model may still be useful in characterizing and comparing 
the damping efficiencies of granular dampers in microgravity 
through the effective frictional force /i M to t ■ 

The effective frictional force of a experimental damper may 
be estimated through the decay of the peak amplitude, as given 
by 



A*M t( 



fc(A - \x n \) 
2n + l - (-l) 2n 



(18) 



where n is the index of the amplitude peak and \x n \ is the 
absolute oscillator displacement for the nth peak. The peak 
number n, used to calculate the effective friction coefficient, 
should be odd to only sample the amplitude during the "col- 
lection" phase of the oscillation and should be as small as pos- 
sible for correct measurement of rapid dampers. The earliest 
value of n which satisfies these requirements is n = 3, at a 
time of t — 1.5/cLi sys tcm- Equation dT8b is used to extract an 
effective friction force for each experimental system and the 
corresponding solutions to Eq. (TTTb are plotted in Figs. 
For the non-optimal dampers the model fits the data well. De- 
viations begin to appear towards the end of the oscillations as 
the "collect and collide" motion begins to break down and the 



granulate spreads uniformly over the box. This lends weight 
to the argument that the "collect and collide" motion of the 
granulate is responsible for the apparent friction-damped be- 
havior. For the optimal damper (see Fig. [7]), the model does 
not fit as well. The discrepancy arises from the oscillators fre- 
quency for this box length being significantly different from 
the predictions of Eq.[T2] Simple spring models will no longer 
work and appear to be constrained to the range (i gi p Ca k S V 4 )- 
A better fit may be obtained by fitting the spring coefficient k, 
however this is unsatisfactory as the results of the model can- 
not be compared between systems. Another deficiency of the 
friction model is that it predicts that the damper will come to 
a complete halt after a finite time. It fails to capture the per- 
sistent small amplitude oscillations (see Fig. |7j). The friction 
deceleration, /.i, still appears to be a useful value for compar- 
ing the damping efficiency of sub-optimal granular dampers. 

The quadratic decay of energy with time in a granular 
system attached to a linear spring has been reported be- 
fore (71 |23l |24j, |47[] . Surprisingly, the same behavior is found 
also for rather different dampers such as thrust-based damp- 
ing l48tl and impact dampers l49l - l52tl . However, this is not a 
general rule and other published results exist (e.g., Ref. (5D) 
where a non-linear decay of the amplitude of the oscillation 
(possibly exponential) is found. This work clarifies that this 
apparent frictional behavior may also arise solely from the 
collisional granular dynamics and does not necessarily arise 
from friction forces within the experimental setup. This is ev- 
ident as the simple model used in the simulations reproduces 
the linear decay of the amplitude. 



VI. CONCLUSIONS 

In this paper, a method for performing controlled experi- 
ments on granular damped oscillators in microgravity is out- 
lined. High-speed video capture and image-processing tech- 
niques are used to reconstruct the motion of the oscillator to 
obtain accurate experimental results. A simple hard sphere 
model and event-driven dynamics are also used to generate 
quantitative results that compare well against the experimen- 
tal values. From the excellent agreement of the simulation and 
experimental frequency, it appears that the damper frequency 
responds like a simple harmonic oscillator to changes in load 
(Eq. ( TT2l ) for short box lengths. This is remarkable given the 
periodic decoupling of the granulate from the spring and box. 
The simulation model scales trivially with the frequency of 
the oscillations as, apart from the negligible initial energy, the 
model has only one time scale. Further research is required on 
experimental systems to determine the frequency dependence 
of granular dampers and generalize the current model to these 
systems. 

The straightforward design of these granular dampers 
yields a remarkably simple expression for the optimal damp- 
ing configuration of the form of Eq. ( TToT ). Simulations at 
the predicted optimal box length damp large amplitude os- 
cillations remarkably well (see Fig. |7ji but are susceptible to 
smaller amplitude disturbances. The final expression for the 
optimal box length is independent of the oscillation frequency, 
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which may be understood through dimensional analysis of the 
model. 

Unlike conventional viscous-damped systems, the granu- 
lar damped system studied here displays a linear decay in the 
amplitude. This behavior is not intuitive and is a feature typi- 
cal of friction-damped oscillators. The simulation results and 
their excellent agreement with experimental results strongly 
suggest that this effect arises solely from the granular dynam- 
ics. The linear decay is a useful property as it implies that 
a granular damper can completely damp oscillations within 
a finite time; however, this is not the case as, at low oscilla- 
tion energies a transition occurs and the damping force is sig- 
nificantly reduced. Further research is required in designing 
dampers with a wider amplitude response by coupling mul- 
tiple dampers with different lengths. The internal geometries 



may also be optimized to eliminate the decoupling of the gran- 
ulate in the midpoint of the stroke to create more effective 
dampers. 
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